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Abstract 

Light quark masses are important fundamental parameters of the Standard Model. The decays ip' — > J /ipiY°{rf) were widely used in determin- 
ing the light quark mass ratio m u /md- However, there is a large discrepancy between the resulting value of 7n u /rrid and the one determined 
from the light pseudoscalar meson masses. Using the technique of non-relativistic effective field theory, we show that intermediate charmed 
meson loops lead to a sizeable contribution to the decays and hence make the ip' — > J/ipir° (77) decays not suitable for a precise extraction of 
the light quark mass ratio. 
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The decays of ip' into J/ipir° and J/tprj were suggested 
to be a reliable source for extracting the light quark mass 
ratio m u /md IHS] (for reviews, see Refs. Orl^l])- The de- 
cay ip' —> J/ipir° violates isospin symmetry. Both the up- 
down quark mass difference and the electromagnetic (em) in- 
teraction can contribute to isospin breaking. However, it was 
shown that the em contribution to the decay ip' — > J/ipn 
is much smaller than the effect of the quark mass differ- 
ence Based on the QCD multipole expansion and the 
axial anomaly, the relation between the quark mass ratio 

= JL (1) 

R m s — m 

where m = (rrid + m u )/2, and the ratio of the decay widths 
of these two decays was worked out up to the next-to-leading 
order in the chiral expansion J8l ISD . At leading order, the rela- 
tion reads Jl^l 
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where q^u,) denotes the momentum of the pion (eta) in the rest 
frame of the ip' and represents SU(3) breaking effects. 
Assuming A^< < 0.4, an upper limit of R was determined 
through Eq. (@]i fioll . It can also be obtained by constructing 
a chiral effective Lagrangian for charmonium states and light 
mesons in a soft-exchange-approximation ll 111 . The ampli- 
tude for the ip' —> J I tpir scales as 



M(ip' -> J/ipir°) ~ (m d -m u )\q w \ 



(3) 



Using the relation between the masses of quarks and 
mesons 1 12, 13], Eq. (f2]i may be rewritten as |0] 
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where and M n r rj \ are the decay constant and mass of 

the pion (eta), respectively. Using Eq. (RJ and the most recent 
measurement of the decay- width ratio [ 14] 



R^o h = (3.88 ± 0.23 ± 0.05)%, 



(5) 



the up-down quark mass ratio is obtained as [ 3 

m 



m d 



0.40 ±0.01. 



(6) 



This value is much smaller than the result obtained from the 



time-honored formula 11211 

m d 



ivi- K+ - M^ + 2M 2 



0.56, (7) 



and it is also smaller than the large N c bound, m u /md ^ 1 /2, 
derived in Refs. H3, E3]. Note that Eq. © is very little 
affected by higher order corrections. It is therefore of fun- 
damental interest to understand theoretically the discrepancy 
between the values of up-down quark mass ratio determined 
from different sources. This Letter is devoted to show that 
the ip' decays into J/ipTr°(ri) are not suitable for extracting 
the quark mass ratio, and hence the seeming discrepancy be- 
tween Eq. (O and Eq. (0 is meaningless. The reason underly- 
ing this statement is that the earlier analysis neglected effects 
from intermediate (virtual) charmed mesons. Those loops 
were shown to be important in some charmonium decays in 



phenomenological models, see, for instance, Refs. IU8I42 1 



As we will show, based on a power counting argument in the 
spirit of heavy quark effective field theory (HQEFT), which is 
supported by an explicit calculation, these contributions over- 
whelm the one directly related to the quark masses. 

To be specific, we calculate the pertinent diagrams for the 
decays ip' — > J/ipir (ri) involving the lowest lying pseu- 
doscalar and vector charmed mesons, see Fig. Q] The cou- 
plings of pion and eta to the charmed mesons follow from 
heavy quark symmetry and chiral symmetry lEH - E^ . In the 
two-component notation of Ref. 112411 . the charmed mesons are 
represented by H a = V a ■ a + P a with V a and P a denoting 
the vector and pseudoscalar charmed mesons, respectively, a 
is the Pauli matrix, and a is the flavor index. The lowest order 



axial coupling Lagrangian is 124tl 



C<t> = ^| Tr [ H tH b a ■ u ba ] 



(8) 
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FIG. 1: The decays ip' — > J/ip-K°(rj) through triangle charmed- 
meson loops. Charmonia, light mesons, pseudoscalar and vector 
charmed mesons, are denoted by double, dashed, thin and solid 
lines,respectively. 



can contribute, and it is the mass differences within the multi- 
plet that generates the isospin or the SU(3) breaking. Similar 
effects have been studied in ao — /o mixing 1E7II28I1 . and the 
isospin breaking hadronic decay of the D* (2317) I29I- 31 1. 
Explicitly, the ip' — > J/ tpir decay amplitude is proportional 
to the difference of the charged and neutral mesons loops 



MW -> J/^tt ) cx e i ^qie^e k J/ ^I c - I n ), 



(10) 



where Ejy^^n denotes the spatial component of the polar- 
ization vector of the J/ip(ip'), I c an d I n are the loop integral 
expressions which will be given below in Eq. ( TT4b for charged 
and neutral charmed mesons. Denoting the expression for the 
strange charmed-meson loop by I s , one obtains the decay am- 
plitude for the ip' —> J I ipr) 



where the axial current is u = -V2d<t>/F + 0(<p 3 ). F is the 
pion decay constant in the chiral limit, and the 3x3 matrix 
<p collects the octet Goldstone bosons. The leading order La- 
grangian for the coupling of the J j ip to the charmed and anti- 
charmed mesons can be constructed considering parity, charge 
parity and spin symmetry. In two-component language, it is 
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with A d B = A{dB) - (dA)B. The charmonium field is 
given by J = ip ■ a + rj c with ip and r/ c annihilating the ip and 
rj c states, and H a = — V a ■ a + P a is the field for anti-charmed 
mesons l26ll . This Lagrangian was first introduced in Ref . ] 25 ] 
in four-component notation with the same coupling <? 2 - Since 
ip' is the first radial excitation of the J/tp, the Lagrangian for 
the ip' coupling to the charmed and anti-charmed mesons has 
the same form as Eq. (0 with the coupling constant 172 re- 
placed by the one for ip', g' 2 . 

Because the ip' and J/ip are SU(3) singlets, it is obvious 
that the decay ip' — > J/ipir violates isospin symmetry, and 
the decay ip 1 — > J/ipr) violates SU(3) flavor symmetry Ir40ll . 
Accordingly, the decay amplitudes reflect the flavor symmetry 
breaking. Here, all the charmed mesons in a flavor multiplet 



k, ^V3 (Ic 



■I n -2I„). (11) 



Before performing the explicit evaluation of the loops it is 
important to first understand the power counting of the sys- 
tem. As was just derived, each vertex in the triangle dia- 
grams is of p-wave character and is thus linear in the mo- 
mentum. Due to parity conservation, one momentum has 
to appear as external parameter (c.f. Eq. (0). Thus the 
loops themselves scale as v 3 /(v 2 ) 2 v 2 = v, where we re- 
placed momentum factors by the dimensionless velocities — 
the proper expansion parameter of HQEFT — and the factors 
denote the non-relativistic integral measure and propagators 
as well as the vertex factors just described, in order. The typ- 
ical heavy meson velocity in the loops may be estimated via 

0.53, where M D is the aver- 

(Mjft + M^)/2. 

The quantities of interest here are differences of loops with 
the remaining terms proportional to m q — this is an energy 
scale of 0(v 2 ). We therefore expect the heavy meson loops to 
scale as m q /v \ q\ which gives some enhancement compared 
to Eq. ®. 

To confirm this power counting estimate and allow for a 
more quantitative statement, we now evaluate the diagrams of 
Fig- 03 explicitly using the non-relativistic technique. Let us 
consider diagram (b) in Fig.[T]as an example of these calcula- 
tions. The decay amplitude in d dimensions is given by 



l(2M D -M^)/M D 
aged charmed-meson mass, and 
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where q is the 7r°(?7) momentum, Ac = Mo* — Mo is 
the mass difference between the vector and the pseudoscalar 
charmed mesons, A(b) = —ax 2 + (c — d)x + d, and 
a = q 2 /4, c = b DD M D + q 2 /2 and d = 2(j, DD .b' DD „. 
Further, fi DD * = M D M D * / (M D + M D -) is the reduced 
mass of the D and D*, b' DD , = M D + M D * — M^> and 
boo = 2Mo — Ej/^ with Ej/^ being the energy of the 
J/ip in the ip' rest-frame. The amplitude in Eq. (TT~2T > has been 

multiplied by a factor of J M^iMj/^MpMp, to account for 
the non-relativistic normalization of the heavy meson fields in 



the Lagrangians given in Eqs. (0 and (|9). The dimensionless 
coupling constants g^oo and g^'OD* are related to the di- 

mensionful ones g 2 and g' 2 via g^oo = fj2 J Mj/^M D and 

9i> r DD» = g'2 yj M^i AI0M0* ■ Note in the last step, we have 
taken d = 4. The integral is finite when evaluated with dimen- 
sional regularization for only a power divergence appears. The 
integral appearing in Eq. (PL?! for c > a > 0, d > 0, which is 
satisfied here, is given by 



r.1 



ix^J A(b) = — I 2a\J c — a + (c — c ) (^/~d — s/c — a^j 



x arctan ■ 



2^/a 2a\[d + (c - c'){y/c^a - Vc 7 ) 



{c-c'Y +4ac' 
2 c + d + 



(c-c') 2 +2a(c'-c) + 4a^c'(c-a) | 3 j~c + \[d 
I 



l + O 



(13) 



We checked that neglecting the a term, which is proportional 
to q 2 , only makes a difference of several percent. Thus, ne- 
glecting the 0(q 2 ) terms, the decay amplitude of any loop 
shown in Fig.[T]scales as 



I = - 




(14) 



where fi(fi') is the reduced mass of the charmed mesons con- 
nected to the J/ij}{^}') in the loop, and 6(6') is the difference 
between the charmed meson threshold and Ej/^(M^i). 

We may now compare the explicit expressions to the power 
counting argument presented above. Since y/2jlb and \J2yJb' 
are approximately the momenta of the charmed mesons in the 
loops, we count them as Mov with Mo and v being the mass 
and velocity of the charmed meson. For the purpose of the 
power counting analysis, one can neglect the difference be- 

/2/j6. Denoting 



tween \/2pb and \f2\jfu . Then one has I 
the charged and neutral charmed meson mass difference by 8, 
we have \i c = fin + 6/2 and 6 C = b n + 25 where the lower 
index corn means charged or neutral. Thus, 



M(i>' -> J/V-tt ) ~ \&\(y/2ii c b e - y/2n n b n ) 
; 2n n + b n /2 



\&\8- 



\Qv\ 



\J2fi n b n 



0{5 2 



(15) 



The mass difference S may be divided into the strong (quark- 
mass difference) and the em contributions as \(rrid — m u ) + 
/3e 2 , see Ref. 13111 . therefore M(ip' — > J/ipn ) scales as 
(nid — m u ) \q v \ jv in line with the estimate given above. 

The validity of Eq. © is based on the assumption that the 
light mesons are produced through soft gluons, and hence at a 



distance much larger than the size of the charmonium, which 
is the basic assumption of the QCD multipole expansion 1132 - 
3411 . Then the matrix element of the soft gluon operator be- 
tween the vacuum and a light meson can be worked out us- 
ing the axial anomaly and chiral symmetry. In the mecha- 
nism considered in this Letter, the light mesons are produced 
through their coupling to the virtual intermediate charmed 
mesons. This kind of mechanism was not included in the 
QCD multipole expansion. These contributions are genuine, 
i.e. there is no underlying double counting. This can be 
seen from the fact that the corresponding integrals are finite 
in dimensional regularization and that the leading terms are 
non-analytic in the quark masses. Comparing Eq. (fT3T > with 
Eq. (O, one sees that the charmed-meson loop effects in the 
amplitude are enhanced by a factor of 1/v ~ 2. Therefore, 
they are more important. 

Assuming the intermediate charmed-meson loop mecha- 
nism saturates the decay widths of the tp' — > J/ipTr°(r]), we 
get 



R* 



'/'/ 



0.14 ±0.09, 



(16) 



where the central value is what we get from a direct calcula- 
tion, and the uncertainty is from neglecting the contribution 

of Eq. © using v - J{2M D - Mi)/M D ~ 0.53 and con- 



tains the one that originates from either using physical masses 
or averaged masses for the field normalizations. This value is 
within 2cr of the experimental ratio. Note that in the ratio all 
the coupling constants g, g 2 and g' 2 disappear. 

We cannot give a prediction for the corresponding decays 
T' — > Yir° (rj) by naively extending the formalism to the bot- 
tom sector. This is because the strong and the em contribu- 
tions to M B o — M B + interfere destructively 13511. and make 
M B o - M B + as small as 0.33 ± 0.06 MeV jlfll, Accord- 



3 



ingly, although the B-meson loop contribution to the decay 
T' — » Try is more important than in the charm sector — v is 
smaller — its contribution to the T' — ^ Yir° is highly sup- 
pressed. On the contrary, in the charmed sector of relevance 
here, the strong and em contributions to M D + — Mpo inter- 
fere constructively J3H1 . and hence enhance the meson loop 
effects. 

Further support of the proposed scheme is provided by ana- 
lyzing the resulting absolute values of the decay widths. Using 
g = 0.6, which is extracted from a tree level calculation of the 
D*+ width, F = 92.4 Me V, and 5 M M = g M ( ,, = 

G, we obtain the absolute values of the decay widths as 

r<y -> J/ipir°) = (3.6 ± 1.9 [5.5 ± 2.9]) x 1CT 4 G 4 keV, 
T(V/ -> J/i)rj) = (2.5 ± 1.3 [4.9 ± 2.6]) x 10~ 3 G 4 kcV, 

(17) 

where the numbers outside and inside the square brackets are 
obtained using the physical and averaged masses for the field 
normalizations, respectively. In order to reproduce the exper- 
imental values r(V>' -> J/ipn") = 0.40 ± 0.03 keV and 
T(ip' -> J/ipn) = 10-0 ± 0.4 keV HI, we need G = 
6.2±0.9[5.5±0.9]andG = 8.4±1.3 [7.1±1.1], respectively. 
These numbers are close to independent model estimates for 
g-ipDD existing in the literature, see e.g. Refs. 1E5II36, 371. 

In this context, we want to stress that extracting the 
quark mass ratio 1/R defined in Eq. (fTJ using T(r)' — > 
7r°7r + 7r~) /T(?/ — > ?y7r + 7r~)has been questioned in Ref. J38ll . 
where meson loops were also shown to be significant. 

In summary, in this Letter, utilizing the technique of non- 
relativistic effective field theory, we show that intermediate 
charmed mesons play an important role in the ip' decays into 
J I ^7T° and J/ipV- They are enhanced by a factor of 1/v com- 
pared with the contribution directly from the quark mass dif- 
ferences. The light quark mass ratio can only be extracted 
from these decays after establishing a complete effective field 
theory up to next-to-leading order with the Goldstone bosons, 
charmonia and charmed mesons as the degrees of freedom in 
the future. What was done in this Letter can be regarded as 
the first step towards that goal. 

We would like to thank B. Kubis for helpful discussion. 
This work is partially supported by the Helmholtz Associa- 
tion through funds provided to the virtual institute "Spin and 
strong QCD" (VH-VI-231) and by the DFG (SFB/TR 16, 
"Subnuclear Structure of Matter"). We also acknowledge the 
support of the European Community-Research Infrastructure 
Integrating Activity "Study of Strongly Interacting Matter" 
(acronym HadronPhysics2, Grant Agreement n. 227431) un- 
der the Seventh Framework Programme of EU. 



[1] B. L. Ioffe, Yad. Fiz. 29, 1611 (1979) [Sov. J. Nucl. Phys. 19, 
827 (1979)]. 

[2] B. L. Ioffe and M. A. Shifman, Phys. Lett. B 95, 99 (1980). 



[3] J. F. Donoghue, Ann. Rev. Nucl. Part. Sci. 39 (1989) 1. 
[4] U.-G. MeiBner, Rept. Prog. Phys. 56, 903 (1993). 
[5] H. Leutwyler, NATO Adv. Study Inst. Ser. B Phys. 363, 149 
(1997). 

[6] J. F. Donoghue and S. F. Tuan, Phys. Lett. B 164, 401 (1985). 
[7] K. Maltman, Phys. Rev. D 44, 751 (1991). 
[8] J. F. Donoghue and D. Wyler, Phys. Rev. D 45, 892 (1992). 
[9] J. F. Donoghue, B. R. Holstein and D. Wyler, Phys. Rev. Lett. 
69 (1992) 3444. 
[10] H. Leutwyler, Phys. Lett. B 378, 313 (1996). 
[11] R. Casalbuoni, A. Deandrea, N. Di Bartolomeo, R. Gatto, 

F. Femglio and G. Nardulli, Phys. Lett. B 309, 163 (1993). 
[12] S. Weinberg, in A Festschrift for I.I. Rabi, ed. by L. Motz, Trans. 

New York Acad. Sci. 38, 185 (1977). 
[13] J. Gasser and H. Leutwyler, Phys. Rept. 87, 77 (1982). 
[14] H. Mendez et al. [CLEO Collaboration], Phys. Rev. D 78, 
011102 (2008). 

[15] J. Z. Bai et al. [BES Collaboration], Phys. Rev. D 70, 012006 
(2004). 

[16] C. Amsler et al. [Particle Data Group], Phys. Lett. B 667, 1 

(2008) and 2009 partial update for the 2010 edition. 
[17] H. Leutwyler, Phys. Lett. B 374, 163 (1996). 

[18] G. Li and Q. Zhao, Phys. Lett. B 670, 55 (2008). 

[19] X. Liu, B. Zhang and X. Q. Li, Phys. Lett. B 675, 441 (2009). 

[20] Y. J. Zhang, G. Li and Q. Zhao, Phys. Rev. Lett. 102, 172001 

(2009) . 

[21] G. Burdman and J. F. Donoghue, Phys. Lett. B 280, 287 (1992). 

[22] M. B. Wise, Phys. Rev. D 45, R2188 (1992). 

[23] T. M. Yan, H. Y. Cheng, C. Y. Cheung, G. L. Lin, Y. C. Lin and 

H. L. Yu, Phys. Rev. D 46, 1148 (1992) [Erratum-ibid. D 55, 

5851 (1997)]. 

[24] J. Hu and T. Mehen, Phys. Rev. D 73, 054003 (2006). 
[25] P. Colangelo, F. De Fazio and T. N. Pham, Phys. Rev. D 69, 
054023 (2004). 

[26] S. Fleming and T. Mehen, Phys. Rev. D 78, 094019 (2008). 
[27] N. N. Achasov, S. A. Devyanin and G. N. Shestakov, Phys. Lett. 

B 88, 367 (1979). 
[28] C. Hanhart, B. Kubis and J. R. Pelaez, Phys. Rev. D 76, 074028 

(2007) . 

[29] A. Faessler, T. Gutsche, V. E. Lyubovitskij and Y. L. Ma, Phys. 

Rev. D 76, 014005 (2007). 
[30] M. F. M. Lutz and M. Soyeur, Nucl. Phys. A 813, 14 (2008). 
[31] F. K. Guo, C. Hanhart, S. Krewald and U.-G. MeiBner, Phys. 

Lett. B 666, 251 (2008). 
[32] K. Gottfried, Phys. Rev. Lett. 40, 598 (1978). 
[33] M. B. Voloshin, Nucl. Phys. B 154, 365 (1979). 
[34] Y. P. Kuang, Front. Phys. China 1, 19 (2006). 
[35] F. K. Guo, C. Hanhart and U.-G. MeiBner, JHEP 0809, 136 

(2008) . 

[36] K. L. Haglin and C. Gale, Phys. Rev. C 63, 065201 (2001). 
[37] A. Deandrea, G. Nardulli and A. D. Polosa, Phys. Rev. D 68, 
034002 (2003). 

[38] B. Borasoy, U.-G. MeiBner and R. NiBler, Phys. Lett. B 643, 41 
(2006). 

[39] There is a discrepancy between the CLEO result ITih given in 
Eq. © and the BES result 0, R„ 0/v = (4.8 ± 0.5)%. If 
we use the branching fractions of the ip' — > J/ipir° and ip' — > 
J/ipV given by the Particle Data Group Ref. [16], which result 
in R n o /v = (4.0 ± 0.3)%, the result m u /m d = 0.39 ± 0.02 
is slightly smaller. 

[40] Here we assume the 77 is in the SU(3) octet, and the effect of the 
77 — 77' mixing is assumed to be small fTTIl . 



4 



